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ABSTRACT 


The  problem  treated  in  this  paper  is  that  of  the 
evaluation  of  pressure  waves  that  can  be  generated  by  clouds 
of  explosive  gas-air  mixtures  in  a  free  atmosphere  which  is 
initially  at  a  uniform  state.  The  particular  subject  of  this  study 
is  concerned  with  the  final  stage  of  the  process  when,  following 
ignition  and  initial  flame  acceleration,  the  flame  acquires  the 
steady  state  of  a  constant  velocity  deflagration.  Under  the 
assumption  of  point,  line  or  plane  symmetrical  geometry,  the 
flow  field  is  then  self-similar.  Theoretical  results  describing 
such  flow  fields  are  presented  and  numerical  solutions, obtained 
for  a  representative  case  corresponding  to  a  hydrocarbon-air 
mixture  at  normal  atmospheric  conditions,  are  given. 
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INTRODUCTION 


When  an  explosive  gas-air  mixture  in  an  unconfined  space  is  ignited  at 
the  center,  either  a  flame  or  a  detonation  can  be,  at  first,  initiated,  depending  on 
the  power  density  of  the  ignition  energy.  Under  normal  circumstances  the  former 
will  then  accelerate  reaching  soon  a  steady  state  of  propagation,  while  the  latter 
will  most  likely  decay  -  a  process  accompanied  by  the  gradual  separation  of  the 
reaction  zone  from  the  shock  front. 

Such  phenomena  have  been  studied  initially  with  reference  to  spherical 

detonations,  by  Manson  Freiwald  and  Ude  ^ \  Zel'dovich  et  al.  ^  and  Leonas^ 

More  recently!  the  knowledge  in  this  field  of  study  has  been  greatly  enhanced  as  a 

consequence  of  the  availability  of  modern  experimental  techniques,  notably  the  use 

of  the  focused  giant-pulse  laser  beam  for  ignition,  as  exemplified  by  the  work  of 
(5) 

Lee;  and  the  application  of  high-frequency  response  pressure  transducers  combined 

with  high  resolution  optical  recording  techniques,  as  illustrated  by  the  studies  of 

(5) 

Soloukhin  and  Oppenheim  '  . 

In  both  cases  of  ignition,  whether  it  results  in  a  flame  or  produces  in¬ 
itially  a  detonation  wave,  the  most  likely  eventual  outcome  of  combustion  in  air  is 
the  establishment  of  a  steady  propagation  state  where  the  flame  moves  at  a  constant 
speed,  acting,  in  effect,  as  a  deflagration  that  leaves  .behind  products  at  rest  and 
generates  ahead  a  pressure  wave  which,  under  proper  geometrical  conditions,  ac¬ 
quires  the  character  of  a  self-similar  blast  wave.  It  is  the  analysis  of  this  case 
that  forms  the  objective  of  this  paper. 

Such  flow  fields  were  studied  by-Taylor  Sedov  ^  and  Shiken 

(Q\ 

and  discussed  in  the  text  of  Courant  and  Friedrichs  '  The  analysis  presented 
here  follows,  in  essence,  the  technique  of  the  above  references;  the  formulation  is, 
however,  modified  to  take  advantage  of  the  simplification!'  that  can  be  obtained  as  a 
consequence  of  our  recent  studies  on  the  fundamental  of  blast  waves!10* 

The  results  are  worked  out  in  detail  to  yield  a  full  description  of  the  pressure  wave 

i 

that  can  be  formed  by  a  cloud  of  a  hydrocarbon-air  mixture.  The  possible  effects  of 
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geometry  are  taken  into  account  by  the  consideration  of  either  point-,  line-,  or 
plane-symmetrical  motion.  In  free  space  the  most  probable  shape  the  wave  may 
take  should  be  encompassed  by  the  point-  and  line-symmetrical  geometry,  so  that 
the  physically  realizable  case  should  he  bounded  by  these  two  solutions. 
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PROBLEM 

The  time-space  wave  diagram  of  the  system  is  presented  in  Fig.  l,the 
coordinates  {  and  f  being  measured  from  the  origin  where  all  wave  fronts  inter¬ 
sect.  The  medium  into  which  the  leading  wave  front  propagates,  state  (p,  as  well 
as  the  products  of  combustion  behind  the  flame,  state  0,  are  considered  to  be  at 
rest.  The  flow  field  to  be  determined  is  self-similar,  so  that  the  trajectories  of 
both  the  deflagration  and  the  front  discontinuity  are  straight  lines,  while  the  gas- 
dynamic  parameters  remain  invariant  along,  similarity  lines  x=  f/Q  *  With  refer¬ 
ence  to  Fig.  1,  the  continuous  self-similar  flow  field  is  thus  associated  with  the 
change  of  state  from  0  to(?).  An  equivalent  flow  field  can  be  generated  by  a  piston 
whose  path  is  also  shown  on  this  diagram.  It  is,  in  fact,  such  a  piston-driven  wave  that 

(Q\ 

formed  the  objective  of  the  first  paper  of  Taylor  '  in  the  subject  of  blast  waves. 

The  problem  is  to  determine  the  space  profiles  of  the  gasdynamic  para¬ 
meters,  i.e.,  the  pressure,^)  ,  density,  ^  ,  temperature,  T ,  and  particle  velocity, U  , 
that  would  be  established  at  any  point  in  space  in  a  combustible  mixture  giving  rise 
to  a  deflagration  front  propagating  into  the  medium  immediately  ahead  of  it  at  a  given 
burning  speed, 5 .  Both  the  deflagrations  associated  with  a  pressure  drop,  as  pre¬ 
scribed  by  the  appropriate  Iiugoniot  curve,  as  well  as  constant  pi'essure  deflagrations 
are  considered  in  the  analysis. 

In  a  self-similar  flow  field  each  particle  undergoes  the  same  thermo¬ 
dynamic  process.  This  is  represented  on  the  pressure-specific  volume  diagram  in 
Fig,  2.  After  compression  by  the  shock  front  from  state (l)  to  state  (f)  represented 
by  a  point  lying  on  the  Rankine-Hugoniot  curve,  RH,  the  particle  undergoes  an  isen- 
tropic  compression  in  the  blast  wave  to  state(p  immediately  ahead  of  the  deflagration 
which,  in  turn,  yields  state  (V)  specified  by  a  point  lying  on  the  appropriate  Hugoniot 
curve,  H.  The  latter  can  be  either  associated  with  a  pressure  drop  terminating  at 
state  (j),  or  it  can  correspond  to  a  constant  pressure  process  ending  at  state  Q), 
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~ •,-*-">-’;^7T^  '■'J^f",'*' •  v?r»7;x.^: 


SOLUTION 


The  self-similar  flow  field  of  a  blast  wave  can  be  described  most  con¬ 


cisely  in  terms  of  the  following  reduced  parameters: 


F*$u  and  Z-(n  «) 


where  JX^ 


d  is  the  local  velocity  of  sound,  while,  with  reference  to  Fig.  1,  ^-cllz/d^  . 

The  reduced  parameters  represent  the  phase  plane  coordinates  of  the 
problem,  so  that  its  solution  is  obtained,  in  essence,  by  determining  the  appropriate 
integral  curve  which  prescribes  /X-jiK n)  ,  subject  to  some  specific  boundary 

conditions.  The  latter  are  given,  on  one  si-.'e,  by  the  Rankine-Hugoniot  relations 
across  the  shock  front  and,  on  the  other,  by  the  jump  conditions  across  the  defla¬ 
gration  front,  yielding  zero  particle  velocity  in  the  burnt  gas  region  immediately 
behind  it.  The  condition  of  F  =  1  corresponds  to  the  state  at  the  piston  face  since 
then,  at  tf=  l£2  ,  one  obtains  XW2^C{  »  i. e.j  the  piston  velocity, XW2t  coincides 
with  the  particle  velocity,  U  ,  In  our  notation,  under  such  circumstances  X-Xp 
while  U-  Up* 

For  a  blast  wave  whose  front  propagates  at  a  constant  velocity,  jJL  -  1. 
If  the  state  of  the  undisturbed  medium  is  uniform,  while  the  substance  is  a  perfect 
gas  with  constant  specific  heats,  the  conservation  equations  for  such  self-similar 
blast  waves  can  be  reduced  to  a  single  differential  equation  governing  the  solution, 
namely: 

2D+J(S-M-F)F  ,2) 

whero  dk)F~  X)4  JZ 

X>=  Z-(t-F)* 

while  J  =  0,  1,  2  for  the  plane-,  line-,  and  point-symmetrical  geometry,  respect¬ 
ively,  and  t  is  the  specific  heat  ratio, 
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Once  the  integral  curve  in  the  phase  plane  is  known,  the  solution  is 
completed  by  the  quadrature 


dJjogx 

dlcgf1 


2) 


D+JZ 


(3) 


while  the  gas  dynamic  parameters  of  the  flow  field  are  determined  from  the 
algebraic  relations 


u-^xF 


(4> 


and,  since  the  flow  is  isentropic, 

-2  -  (A}2-  (JL\  W~/-P  I  ~  X 

%  k)  'If,)  -($  -  *  z, 


2  JZ 


At  F-  1  Eqs.  (2)  and  (3)  become  simply: 


(5) 


(6) 

(?) 


dim?  2 

ClltqF  ~  j+i 

Ak.  *  _  _j- 

ddcgF  Ft 

representing  then  also  the  governing  differential  equations  for  the  plane-symmetri 
cal  case  of  J  =  0.  The  solution  in  this  case  is  thus  obtained  immediately  by  inte¬ 
gration,  giving 


Z=ZPF2 


and 


X^XfF 


-/ 


(8) 

(9) 


Zp-ZgFxl  and  . 
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«*I4. 


The  solution  to  the  problem  is  given  in  terms  of  a  family  of  integral  curves 
on  tire  phase  plane  for  which  ^  and  pare  the  coordinates.  Each  curve  corresponds 
to  a  fixed  value  of  Zy  and  Xp  .  They  are  shown  for  the  case  of  /  =  1.3,  on 
Fig.  3,  labeled  in  terms  of  jfZZp  ,  for  spherical,  cylindrical  and  plane  flow 
fields  corresponding,  respectively,  to  J  -  2,  1  and  0. 

The  curves  were  obtained  by  the  numerical  integration  of  Eq.  (2)  or, 
for  J  =  0,  directly  from  Eq.  (8),  starting  from  the  conditions  at  the  piston  face, 
that  is  F  =  1  and  Z  =  Zp  =  t/$  .  They  were  terminated  at  the  shock  front  whose 
phase  plane  coordinates  are  given  by 

. 

and  (10) 

2  (m2l~Fr~yK'?-+yj 

where  y5  ^‘Y^t2-  YM2  •  Eliminating  the  latter  from  £qs.  (10)  and  (11)  one  ob- 
tains  for  the  locus  of  end  points 


representing  the  equation  of  the  Rankine-Hugoniot  curve  in  the  phase  plane.  The 
corresponding  initial  condition  is  that  of  a  quiescent  atmosphere  ahead  of  the  wave, 


that  is 


ff=o 

#.-y 


while  the  change  of  the  gasdynamic  parameters  across  the  fi’ont  is  specified  in  terms 
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of  the  following  shock  relations: 


h-H 

p,~m 


(fr-v) 


-  2M/J1_ 

mui 


«m  <■ 


>) 


(13) 


For  weak  waves  .the  Mach  number  of  the  front  approaches  unity  and 
there  arises  a  basic  difficulty  in  the  integration  of  the  differential  equation  govern¬ 
ing  the  structure  of  the  blast  wave.  This  is  due  to  the  particular  property  of  the 
singularity  of  the  point  F=  0,  2!  ~  1  which  lies  at  the  intersection  of  the  Rankine- 
Hugoniot  curve  with  the  p  =  0  axis  .corresponding  to  the  condition  of  M  =  1  at  the 
front  The  problem  is  caused  by  the  fact  that  the  integral  curves  representing 
the  blast  wave  structure  intersect  the  Rankine-IIugoniot  curve  at  a  positive  slope 
tliat  becomes  larger  for  weaker  waves,  while  they  have  to  approach  the  singularity 
at  F=  0  with  an  infinitely  negative  slope.  The  consequent  rapid  change  in  the  deri¬ 
vatives  produces  queer  numerical  results  which  led  Taylor  ^  initially  to  the  belief 
that  there  may  be  a  lower  bound  to  the  solution  corresponding  to  the  condition  that  X 

i 

cannot  be  smaller  than  a  certain  value,  somewhat  between  0,4  and  0.5. 

Actually,  as  it  has  been  indeed  observed  by  Taylor  when  the  front 
Mach  number  approaches  unity,  the  problem  is  governed  by  the  acoustic  wave  equa¬ 
tion  and,  as  it  will  be  demonstrated  here,  there  exists  a  continuous  transition  be¬ 
tween  the  acoustic  and  blast  wave  domains,  so  that,  in  fact,  the  weak  limit  does  not 
exist  and  X  can  be  small  as  one  pleases. 

The  acoustic  wave  is  described  in  terms  of  the  velocity  potential 
which  satisfies  the  equation 

Hi'  ~  ri  m  arJ 


(14) 


and  is  related  to  the  particle  velocity  and  pressure,  respectively,  as  follows 

-34 

tf'  ‘fist 


(15) 


Introducing  the  similarity  variable  ,  and  assuming  for  the 

potential  the  functional  foim  (f>(  ft  h )  -  ,  Eq.  (14),  is  transformed  into  an 

ordinary  differential  equation 


x(i-x2)-f"  jf'  =  O 


(16) 


while  Eq.  (15)  yields 

“  Y! •  i-[f« -*?“>]'  rlf.  -  9>] 

Integrating  Eq.  (16),  subject  to  the  boundary  conditions  at  the  piston 

face  where 


~  ~Q,  $ 

(18) 

and  at  the  front  where 

fsd><  at 

''io 

11 

(19) 

one  obtains 

'-Xz  j 7/2 
-Xf) 

UtxL 

m  =  ax  -W  (i 

(20) 
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and 


The  corresponding  velocity  and  pressure  profiles  can  be  then  express¬ 


ed  as 


(22) 


for  j=l 


for  j=2 


(23) 


The  above  set  of  equations  provides  a  complete  description  of  the 
acoustic  wave  in  terms  of  a  single  variable  y  ,  the  self-similarity  coordinate,  and  the  para 
meter  equivalent  to  ,  as  it  can  be  verified  directly  from  the  defini¬ 

tions  of  %  and  Zp  by  noting  that,  for  the  acoustic  wave,  y  =  1  while  Op^Ql  . 

For  the  plane-symmcti'ical  case  (j  =  0)  the  acoustic  wave  is  identical 
to  the  blast  wave  corresponding  to  y  =  1,  one  being  just  a  simple  limit  of  ;he  other 
without  any  anomaly. 

The  flame  is  considered  to  act  as  a  deflagration,  as  depicted  on  Fig,  2. 

From  the  analysis  of  the  thermodynamic  processes  represented  on  this  diagram, 
taking  properly  into  account  the  dependence  of  the  heat  of  reaction  on  the  tempera¬ 
ture  immediately  ahead  of  the  flame, one  obtains  the  relation 


vF = f = x'ilik  -  f£o*  &?)]-  (fyr  }p§-i)  w 


where  is  the  molecular  weight  ratio. 
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In  the  course  of  computations  the  parameters  Zt /*■  and  X/  are 
evaluated  point  by  point  starting  from  the  piston  face.  The  corresponding  value  of 
))p  can  be  therefore  determined  from  Eq.  (24)  and  the  position  of  the  flame  front  is 
thus  established  when  it  attains  the  appropriate  magnitude  for  a  given  combustible 
mixture. 

The  two  terms  in  the  square  brackets  represent,  respectively,  the  two 
components  of  the  expression 


Zi  (H-p) 

while 


(25) 


(26) 


For  a  constant  pressure  deflagration,  =  and,  consequently,  the 
second  term  in  the  brackets  vanishes.  If,  moroover,  the  influence  of  heating  in  the 

t'Zi-lQi 

W  •"#  r  «...  .'ikJ  f 

form 


blast  wave  is  negligible,  X* -/^iz £  {  and  Eq.  (24)  is  reduced  to  the  asymptotic 

( Qt!  ~ 


(27) 


The  value  of  the  parameter  /-=£•  at  the  flame  front  is,  under  such  circumstances,  es¬ 
tablished  solely  on  the  basis  of  the  value  of  Vp  as  specified  above. 

Finally  it  should  be  noted  that  to  each  point  of  the  integi-al  curve  there 
corresponds  a  certain  value  of  the  burning  speed  given  by  the  relation: 


5=  Xi^z(l-Ft-) 


(28) 


RESULTS 


For  numerical  computations, a  typical  case  of  a  hydrocarbon-air  mix¬ 
ture  initially  at  N.  T.P.  conditions  was  used.  Its  thermodynamic  properties  were 
expressed  in  terms  of  the  following  parameters 

Vp  =  7  Cl  i  =  345  /sec  ^  ft1  ~  1  atm 

=  =  1.2  %- i 

Figure  3  represents  integral  curves  on  the  phase  plane  together  with 
the  corresponding  plot  of  .  Each  is  labeled  with  the  appropriate  value  of 

the  particular  /  Zp  .  Points  corresponding  to  the  flame  front  are  indicated 
by  triangles.  As  pointed  out  before,  the  transition  between  the  blast  wave  solution 

r2 

-  10  the  two  are  exactly  coincident. 

•Y  2 

For  s  >  10  the  acoustic  solutions  become  quite  accurate,  while  the  numerical 
computations  for  the  blast  wave  get  to  be  very  lengthy  and,  hence,  suffer  from  loss 
in  accuracy. 

To  give  a  complete  physical  description  of  the  flow  fields  thus  evaluated, 
the  corresponding  particle  velocity,  temperature  and  pressure  profiles  for  the  spher¬ 
ical  case  of  j  =  2  are  presented  in  Figs.  4,  5  and  6,  respectively. 

Figure  7  gives  the  burning  speed,  S  ,  evaluated  from  Eq.  (28),  as  a 
function  of  Vp  determined  by  the  use  of  Eq.  (24)  for  the  spherically  symmetric 
waves,  with  £  as  the  parameter.  The  graph  is  quite  general  in  that  it  is  applicable 
to  any  value  of  the  heat  of  reaction  that  may  be  ascribed  to  the  combustible  mixture. 
As  on  the  other  diagrams,  points  corresponding  to  the  given  substance  for  which 
^p=  7  are  denoted  by  triangles.  The  burning  soeed,  S  ,  has  an  upper  bound  given  by 
the  CJ  deflagration  limit.  It  is  specified  by  the  condition  ^=#4  and  is  repre¬ 
sented  on  the  diagram  by  the  chain-dotted  line. 
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The  most  significant  physical  parameters  for  all  the  geometries  are 
shown  in  Fig.  8,  plotted  as  a  function  of  %  ,  the  parameter  identifying  specific 
solutions.  This  includes  the  burning  speed,  S  ,  and  the  flame  position,  ,  both 
referring  to  a  typical  hydrocarbon-air  mixture  that  corresponds  to  \)p  =  7,  as 
well  as  the  front  parameter  y  =  J/Ma .  The  latter  exhibits  most  distinctly  the  effect 
of  the  geometry  on  the  strength  of  the  wave  generated  by  the  flame  —  the  Mach 
number  increasing  considerably  from  a  point-  to  a  plane-symmetrical  case  for  the 
same  value  of  S  ,  while  the  corresponding  change  in  Xj  is  quite  insignificant. 

Finally  Fig.  9  displays  significant  pressure  levels  as  a  function  of  the 
burning  speed  S  attained  in  a  combustible  medium  corresponding  to  \^:  =  7.  For 
the  same  burning  speeds,  the  differences  between  the  three  geometries  are  quite 
significant,  the  plateau  of  for  the  plane-symmetrical  wave  being  appre¬ 

ciably  higher  than  the  pressure  levels  in  the  other  two  cases.  In  practice  one  can 
expect  such  high  pressures  to  be  attained  whenever  the  flame  front  is  well  confined, 
as  in  corridors  or  tunnels. 
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